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ABSTRACT 



We collect and combine all published data on the vertical atmospheric muon flux and the 
muon charge ratio for muon momenta above 10 GeV. At sea level the world average of the 
momentum spectra agrees with the flux calculated by E.V. Bugaev et al. within 15 %. The 
observed shape of the differential flux versus momentum is slightly flatter than predicted in 
this calculation. The experimental accuracy varies from 7% at 10 GeV to 17% at 1 TeV. The 
ratio of fluxes of positive to negative muons is found to be constant, at a value of 1.268, with 
relative uncertainties increasing from approximately 1% at low momenta to about 6% at 300 
GeV. 
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1. Introduction 



We collect measured atmospheric muon flux data and charge ratios as a function of momen- 
tum and compute world averages. Only measurements at sea level or low altitudes and for 
(near) vertical incidence are taken into account, since several data sets are available for these 
experimental conditions. Here we consider only data with muon momenta above 10 GeV. At 
lower momenta geomagnetic effects and solar influences play a significant role and make the 
interpretation of the data more difficult. A recent compilation of charge ratio data at low muon 
momenta can be found in reference 

A precise knowledge of the muon spectrum and charge ratio allows to constrain the primary 
flux and the models of atmospheric showers so that also the atmospheric neutrino fluxes can 
be calculated with a good precision. This is a very important issue, since the Superkamiokande 
experiment^] and others have seen indications for a disappearance of atmospheric muon neu- 
trinos. So far this interpretation is based on the angular distribution and on the ratio of muon 
neutrino to electron neutrino fluxes. It is very important to compare also directly the measured 
and calculated absolute muon neutrino fluxes; until now this was prevented by the large model 
uncertainties. 

2. Effects relevant for spectrum and charge ratio 

The following effects might influence the measurements of the muon flux and the charge ratio. 
It is possible that the published data need to be corrected accordingly in order to arrive at a 
meaningful comparison between the various measurements. 

• Geomagnetic effects 

For near vertical incidence the geomagnetic cutoff for primary protons is below 10 GeV for 
all latitudes at which the cosmic ray measurements were made0 (exceptions are discussed 
below). Geomagnetic effects can therefore be neglected. 

• Solar modulation 

Using the parameterization given in reference {§] we estimate that the primary proton 
flux at 50 GeV (100 GeV) decreases by 3% (1.6 %) at maximum solar activity compared 
to the minimum. The mean primary proton momentum resulting in 10 GeV muons at 
sea level exceeds 100 GeV. Using the air shower program CORSIKA ||, we found that 
about 80% of those protons have a momentum larger than 50 GeV. This results in an 
uncertainty of ±1% for the muon flux at a momentum of 10 GeV. Similarly, one can 
estimate a flux uncertainty of ±0.5% at 20 GeV and less at higher momenta. At 10 GeV 
the charge ratio is expected to change by about ±0.2%. At higher momenta the effect is 
even smaller. We do not correct the data for time dependent solar effects. 

• Altitude dependence 

Not all experiments measure at sea level. In order to investigate the dependence of flux 
and charge ratio on the altitude we used the air shower simulation program CORSIKA 
and also apply the empirical formula found by De Pascale et al.H. 

For muon momenta above 10 GeV and altitudes less than about 1000 m the vertical muon 
flux can be parameterized by 



e h/L ± 0.003 



(1) 



$(h = 0) 
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where h = altitude, L = 4900 m + 750 m and p = muon momentum. 

The form of the parameterization is similar to the one used in||. The uncertainty of 
±0.003 reflects the quality of the parameterization and the comparison to the mea- 
surements. Example: For h = 1000 m and p = 10 GeV we obtain the flux $(/t) = 
1.08 • = 0). Note: Caprice data0 disagree with both @ and CORSIKA for higher 
momenta; here they have not been taken into account. 

The charge ratio is not affected, it changes by less than 0.005 for h < 1000 m and p > 
10 GeV. 

We do correct all published fluxes using formula (1). 
Zenith angle dependence 

The muon data are normally collected within a certain cone around the vertical direction, 
including zenith angles up to Q max . With help of CORSIKA we find that the zenith angle 
dependence can be parameterized in the form 

d$ 

l + a(p) • (l-cos0) (2) 



dcos^ 



with a momentum dependent coefficient a(p). Accordingly we estimate the following flux 
reduction factors 

d$ d$ 

9(0) = ~, -Ml- -(0) (3) 

d cos 9 d cos 9 



p/GeV 


a 


0(5°) 


0(10°) 


<?(20°) 


10 


-1.50 


0.994 


0.978 


0.910 


30 


-1.28 


0.995 


0.981 


0.925 


100 


-0.94 


0.996 


0.986 


0.944 


300 


-0.61 


0.998 


0.991 


0.963 


1000 


-0.22 


0.999 


0.997 


0.987 



Note: the entries are differential values, they have not been integrated over 9. 

Since not all experiments quote the range of accepted zenith angles, we cannot correct 
for this effect. We have to keep in mind that this might cause a bias, especially at low 
momenta. 

Atmospheric pressure/temperature profile 

Previous calculations [|8[ and measurements indicate that the relative muon flux vari- 
ation A$ at ground level is related to the temperature-distribution in the atmosphere 
via 

$ 1 eff 

$ is the integral flux above a certain muon momentum threshold pth- T e fj is the absolute 
effective temperature of the higher atmosphere, a is the temperature coefficient, which 
is a function of zenith angle and muon energy. For zenith angles 9 ~ OH, M: 



a 



70 GeV 
1 + 

Pth 



(5) 
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Simulations using CORSIKA arrive at similar conclusions. Example: At a threshold of 
70 GeV the formula yields a = 0.5. Since the atmospheric temperature, with a typical 
value of 220 K, varies over the year by up to ±5 K, this implies a flux change of ±1%. 

For muon momenta above 10 GeV the pressure at ground level is not expected to show a 
significant correlation with the fluxpl]. 

Unfortunately, most experiments do not report the atmospheric temperature, nor do they 
correct for this effect. It is even not clear, how to define the reference value. Therefore, 
we cannot correct for atmospheric effects. 

Unfolding of the momentum spectrum 

The measured muon spectrum agrees with the true spectrum only if the momentum 
resolution is small compared to the momenta being investigated. Otherwise, the steepness 
of the spectrum, which falls off approximately according to 

d$ o 

(6) 

leads to an asymmetric distortion, an enhancement of the measured flux at high momenta. 
Thus, the measured spectrum needs to be unfolded for experimental resolution effects. 
In the most simple approach - assuming the spectrum is roughly known - this can be 
achieved by a simple correction factor, which has been calculated in JTIJ. The authors 
assume the spectrum (^) and a Gaussian error distribution in the variable 1/p with width 
o~i/ p . Often the experimental resolution is given in terms of the 'Maximum Detectable 
Momentum' Pmdm, defined as the momentum value for which the integral over the Gauss 
distribution becomes 1/2: 



V2 



(Jl/p 



Thus, 



- with E(x) = 4= / e"* 2 dt (7) 
2 V 71 " Jo 



0.6745 cti /p (8) 



Pmdm 

The ratio of the measured and true spectra is then given by 

R (^£) = £(0.4769 **™L) + 1.1829 exp f -0.2275 %M)^ ( 9 ) 

P P Pmdm \ P J 

The measured flux must be multiplied by 1/R to correct for the experimental resolution. 
Figure [l| shows the dependence of R on 1/p. For p < 0.3 Pmdm the correction amounts 
to less than 1% and can be neglected. For higher momenta the correction rises strongly 
and must be taken into account. 

We have assumed that the experimenters have corrected their data for momentum reso- 
lution effects or that they can be neglected. However, several papers are not very clear 
on this point. Therefore, some published spectra might be biased towards too high flux 
values at large momenta. 
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Figure 1: Unfolding correction factor 



3. Experimental Data 

Only published results are taken into account. In appendix A we summarize the characteristics 
of all relevant experiments/publications, in chronological order. The spectrum and charge ratio 
data used in this compilation are listed explicitly in appendices B and C. 

4. Absolute muon flux 

There are two aspects to the measurement of the absolute muon flux, namely the shape of 
the spectrum as a function of energy and the absolute normalization. Some experiments only 
measure the relative muon flux as a function of momentum. Therefore we will analyze the data 
in two steps. First we check the spectral shape, leaving the normalization as a free parameter. 
Secondly we determine the absolute normalization of the spectrum. 

4.1 The shape of the muon spectrum 

A whole range of experiments are performed to measure the muon flux, the measurements used 
are listed in appendix B. We have corrected the datasets for altitude, which is a small correction 
in most cases. In order to be able to compare the datasets, we fit each set to a reference shape, 
using the data with momenta above 10 GeV. In this fit, and in the following, we assume that 
the measurement performed in each momentum bin is independent of the other momentum 
bins. The reference shape is taken from the theoretical calculation by Bugaev et al Jll| , leaving 



the normalization as a free parameter. In general this shape provides a good description of the 
datasets, as can be seen below. The results of our fit are listed in table ED. In this table, we 



separated the data from Rastin[p25[ into two sets, as different normalizations are used in their 
paper. Next to fitting the normalization, we also calculated the normalization by comparing 
the integrated flux above 10 GeV/c (or the lower cutoff of the experiment whichever is higher) 
to an integrated flux calculation using the reference shape. The data published in references 



13 1 , [16] and [15| do not allow for this normalization method. As these papers are normalized 
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Data set 



x 2 /ndf 



Normalization 
from fit 



Normalization 
from integration 



Caro 1950 
Owen 1955 pfl 



Pine 1959 111 



Pak 1961 p| 
Holmes 1961 M 



Hayman 1962 |JT7| 
Aurela 1967 



Appleton 1971 @ 



Allkofer 1971 |0 



Bateman 1971 M 



Nandi 1972 H 



Ayre 1975 f2| 
Green 1979 §| 
Rastin 1984 
Rastin 1984 J25J 
De Pascale 1993 



10-25 GeV) 
>25 GeV) 



Tsuji 1998 26 



Kremer 1994 data [f27j 
Kremer 1997 data [0 



2.6/4 
0.5/2 
4/11 

4/6 
43/12 
13/14 
0.7/2 
38/23 
116/8 

8/8 
60/14 
348/44 
2.3/4 
0.2/5 
24/29 

7/5 
16/13 

10/6 

13/6 



0.65 ± 
0.819 ± 

0.76 ± 

0.75 ± 
0.807 ± 
0.735 ± 

0.81 ± 
0.370 ± 
1.058 ± 
0.871 ± 
0.998 ± 
0.980 ± 

0.98 ± 
0.995 ± 
0.960 ± 
0.798 ± 
0.961 ± 
0.822 ± 
0.831 ± 



0.03 

0.013 

0.03 

0.03 

0.016 

0.007 

0.03 

0.003 

0.006 

0.008 

0.008 

0.002 

0.02 

0.003 

0.005 

0.016 

0.014 

0.009 

0.008 



0.74 ±0.09 
0.829 
0.76 
0.76 
0.829 
0.746 ± 0.008 

0.79 ±0.03 
0.366 ±0.003 
1.01 ±0.01 
0.83 ±0.03 
1.001 ±0.008 
0.95 ±0.02 
0.98 ±0.02 
0.977 ±0.002 
0.951 ±0.005 

0.80 ±0.03 
0.972 ±0.014 
0.818 ±0.007 
0.821 ±0.007 



Table 1: Normalization of datasets with respect to the Bugaev calculation 



to Rossi we recalculate this normalization point with the reference shape. The data of JH 
are normalized to the differential flux at 10 GeV/c. The results of this calculation are shown 
in the last column of table [I]. In general, both normalizations are in good agreement. 

The high energy part of the Rastin data is shown in figure |2|. Here and in the following we 
present all spectra weighted with p 3 , a common practice to compensate for the steep fall-off 
with momentum. Figure § nicely shows that the reference shape fits the data rather well, which 
justifies the use of the Bugaev curve as a reference. However, the data has the tendency to be 
slightly higher than the normalized curve at the higher momentum values. 
The x 2 °f the fit as listed in this table made us re-check five datasets; the first four are shown 



in figure |3|. The data of Holmes |16j clearly show that a simple re-normalization will not work. 



The data points do not follow the reference shape, especially at higher momenta. In their paper 
Holmes et al. apply additional corrections to the highest two data points, indicating that these 
are close to the MDM of the detector. Unfortunately, the value of this maximal momentum is 
not mentioned. 

The Appleton data[l9| are scattered a lot around the curve. With the value of x 2 /NDF being 



only slightly less than 2, this plot suggests that the errors could be underestimated. 
The Allkofer data[fj0] have a completely different shape. The data rise faster than the reference 
shape and plateau at a lower value. This plateau also seems to be wider than suggested by the 
reference distribution. 



The Nandi data|29[ rise to a significantly higher value than predicted by the reference shape. 
This and the low value of the third data point create the large x 2 ■ 

The Ayre data p3) (left side of fig. f|) start off below the curve and continue to rise longer than 
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Figure 2: Muon flux data by Rastin et al.|25[ in comparison to the reference spectrum from 
Bugaev et al. ||llj| , after normalization 
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Figure 3: Muon flux data by Holmes et al.|16|, Appleton et al.[]nj, Allkofer et al.[2Tl] and Nandi 



et al.[29] in comparison to the reference spectrum from Bugaev et al.[ll], after normalization 
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Figure 4: Muon flux data by Ayre et al. ||23|| in comparison to the reference spectrum from 
Bugaev et al.|lT|) after normalization. Left: data as published. Right: spectrum after momen- 
tum scaling. 
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expected. Therefore, the peak is at a higher value, but the drop-off rate seems to be similar as 
predicted on a log(p) scale. The curve suggests that the momentum could be over-estimated. 
A best fit of the momentum scale leads to a scaling of the momenta by a factor of 0.825, see 
right side of fig. |j. The x 2 /NDF improves from 348/44 to 136/44. It naturally changes the 
normalization. The muon spectrum closely follows a p~ 3 -dependence, thus the normalization is 
changed to about 56 % of the original. This is in fact what we observe. Even using the modified 
momentum, for which we cannot find a justification, the Ayre data do not fit the curve very 
well. 

If we ignore these five datasets for the moment, we can compare the remaining data to the 
reference curve. We do this by applying the normalization calculated as outlined before, and 
listed in table [l] as 'normalization from integration'. The result is shown in figure [|. The top 
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Figure 5: 'Good' data in comparison to reference spectrum. 



part of this plot shows a direct comparison of all the datasets to the theory, while the bottom 
part shows the relative difference between these measurements and the description. The larger 
differences at higher momenta are mainly due to the data by Rastin et al. p5| . If we include 
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Figure 6: All data in comparison to reference spectrum. 



the remaining five datasets the difference between the shape of the data and the reference shape 
increases, especially at higher momenta. This can be seen in figure |6|. 

The relative differences are plotted in 8 bins per decade, equidistant in log(p). This choice 
represents a good compromise taking into account experimental uncertainties and the rate of 
change of the p 3 -weighted spectrum with momentum. The y 2 /NDF values and the average 
flux values in the units m -2 sr _1 s _1 GeV -1 are shown for each bin in table 0. The five datasets 
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Table 2: Average flux in m 2 sr 1 s 1 GeV 1 

discussed above have a large impact on the \ 2 of the relative difference in these bins. Therefore, 
we will exclude them when adjusting the shape according to the measurements. We fit a third 
degree polynomial to the logarithm of the flux as a function of the logarithm of momentum. 
We parameterize this function as follows: 

H(y) = H x • (y 3 /2 - 5y 2 /2 + 3y) 

+ H 2 ■ (-2y 3 /3 + 3y 2 - 10y/3 + 1) 

+ H 3 ■ (2/ 3 /6 - y 2 /2 + y/3) 

+ S2-(y 3 /3-2y 2 + lly/3-2) 
y= 10 log(p/GeV) 
F(p) = 10 H{y) m.- 2 si' 1 s- 1 GeV- 1 

This parameterization is similar to the one used by pT| , however the fit variables are chosen 



such that they have a simple interpretation: Hi, H2, and H3 represent the logarithm of the 
differential flux at 10, 100 and 1000 GeV, S2 represents the exponent of the differential flux at 
100 GeV. The \ 2 /NDF of this fit is 8/12, the correlation matrix is shown in appendix D. The 
fitted parameters are listed below, as well as the equivalent values from our reference shape. 
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-5.76 ±0.03 


-2.10 ±0.03 




P,(GeV) 

Figure 7: The ratio between the normalized flux data and the fitted function. The band 
indicates the uncertainty of the fit. The dashed curve represents the normalized reference 
shape |TT . 



Figure ^ contains a la error band, which can be approximated by 

P 



0.003 ±0.00015 • 



GeV 



(11) 



It represents the size of the combined relative experimental uncertainties as a function of mo- 
mentum bin. 

The shape uncertainty at reference momenta of 10, 100 and 1000 GeV are 0.5%, 1.8% and 
15% respectively. Above 200 GeV the uncertainty rapidly increases, indicating that more 
measurements are needed at these momenta. 

4.2 The absolute normalization of the muon spectrum 

In section 4.1 we fitted a functional shape to all renormalized datasets. The renormalization 
was such that the integral flux above 10 GeV corresponds to the calculation by Bugaev. We 
will now fit the functional shape obtained in 4.1 to the datasets of those experiments providing 
an absolute flux measurement, while leaving the normalization as a free parameter. Therefore 
we fit the function 



F(p) =C-10 H{y) 



(12) 
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Data set 


X 7NDF 


Fitted Normalization C 


Allkofer 1971 p0| 
Ayre 1975 23 j 
Bateman 1971 [21] 
Green 1979 24\ 
Tsuji 1998 [26| 
De Pascale 1993 § 
Kremer 1994 data [27 
Kremer 1997 data H 


117/8 
286/44 

8/8 
2.4/4 
16/13 

7/5 

9/6 

11/6 


1.043 ±0.006 
0.964 ±0.002 
0.860 ±0.008 
0.967 ±0.022 
0.948 ±0.014 
0.787 ±0.015 
0.811 ±0.009 
0.820 ±0.008 



Table 3: Normalization factor C with respect to the integrated flux as calculated by Bugaev 

The result is shown in table ^. We again ignore the datasets with a very high x 2 /NDF (Allkofer, 
Ayre). The three remaining data sets with the largest normalization factors (Bateman, Green 
and Tsuji) are measurements performed with solid iron magnet spectrometers, whereas the 
other three (Kremer 1994 and 1997 and De Pascale) use the same superconducting magnet. 
We will first average the normalizations performed by the same collaboration (Bateman and 
Green, Kremer), and afterwards calculate the normalization measurements performed by the 
solid iron magnet spectrometers and the superconducting magnet spectrometers, which gives 
the following results: 



Solid iron magnets 
Superconducting magnet 



0.937 ±0.012 
0.811 ±0.007 



These values are clearly not in agreement. We have no explanation for this observation. We 
will simply take the average of these two values to be our normalization and half the difference 
to be the uncertainty. We arrive to our final value of 0.874 ± 0.063, thus a normalization with 
a relative uncertainty of 7 %. 

4.3 The muon spectrum 

In the preceding sections we have parameterized the muon spectrum at sea-level. We summarize 
the parameters obtained: 

C H\ H2 H% S2 

0.874 ±0.063 0.135 ±0.002 -2.529 ± 0.004 -5.76 ±0.03 -2.10 ±0.03 

Our description, as well as the calculation from Bugaev, and the measurements used for the nor- 
malization are shown in figure |8|. The error is given be the estimated normalization uncertainty 
of 7% and the shape error in (|TT|) , added in quadrature. 



13 




P,(GeV) 



Figure 8: The result on the muon flux. The dotted lines show the 1 sigma error band, whereas 
the dashed curve is the description by Bugaev. The points are the data used in the normalization 
procedure. The open points stand for experiments using a superconducting magnet, the black 
points indicate conventional magnets. 
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If we compare our description of the differential flux to the theoretical description of Bugaev 



lTJ we get the following : 



Momentum: 10 GeV 100 GeV 1000 GeV 

Our Result 1.19 ±0.08 (2.59 ±0.19) • 10~ 3 (1.52 ± 0.26) • 10~ 6 m~ 2 sr" 1 s" 1 GeV" 1 
Bugaev et al. 1.34 2.89 • 10~ 3 1.39 • 10" 6 m~ 2 sr" 1 s" 1 GeV" 1 

At 10 GeV the measured flux is 89% of the calculation of Bugaev. However, the measured 
shape is slightly less steep, and at 1 TeV we arrive to a value which is close to the predicted 
one. 

5. Charge Ratio 



The charge ratio _R M is defined as the ratio of vertical fluxes for positive and negative muons at 
sea level. 

The measured charge ratios together with the published uncertainties are listed in appendix C. 
Figure ^ shows all values as a function of momentum. 




Figure 9: Measured charge ratios. 



In order to study the momentum dependence we have grouped all 15 data sets into momen- 
tum bins chosen to be equidistant in log p. We have combined the different measurements by 
assuming that they are uncorrelated. The bin size is relatively large, since a strong momentum 
dependence is not expected. The result is shown in figure [Hj The two data points around 
500 GeV and the single measurement above 1 TeV have huge uncertainties (~ 20%) and are 
therefore not included in the figure. 

For all eight momentum bins the \ 2 values are good or at least acceptable; this implies the 
various experimental data agree among each other. 
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Figure 10: Average charge ratio values. 



The third momentum yields the highest charge ratio. However, this result can not be attributed 
to a single 'outlier'. 

To see if there is a momentum dependence we have performed the following fits to the charge 
ratio values shown in figure [H| as a function of logp: 

a) /(logp) = = const 

This gives a good fit with x 2 /NDF = 158/142. The resulting charge ratio of 

Rl = 1.270 ± 0.003 (13) 

is displayed in figure |H] as horizontal line. 

b) /(logp) = Rl + Si ■ log(p/GeV) 

Naturally this fit is satisfactory, too. The slope comes out as 

Si = 0.006 ± 0.011 (14) 

which is compatible with zero. 

Therefore, the measured charge ratios are consistent with the hypothesis of being momentum 
independent in the range 10 GeV < p < 300 GeV. 

We have looked at the data in more detail and tried to answer the following questions: 

i) Do the different experiments agree with each other ? 

ii) Is the 'peak' at about 30 GeV significant ? 

The previous statistical analyses and figure ||] seem to imply the answer 'yes' to the first question. 
However, when separately plotting the two (by far) most precise data sets (Baxendale 1975[PO 



and Rastin 19840), one finds the discrepancy displayed in figure |TT[ Averaged over all 



momenta the mean values 

R\{Baxendale) = 1.2799 ± 0.0042 (15) 
R° (Rastin) = 1.2511 ± 0.0058 (16) 
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Figure 11: Comparison of charge ratios as measured by Baxendale et al.|30| and Rastin et 
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20.0 


1.264 ± 0.009 


1.4-1.6 


31.6 


1.283 ± 0.011 


1.6-1.8 


50.1 


1.265 ± 0.014 


1.8-2.0 


79.4 


1.252 ± 0.017 


2.0-2.2 


126 


1.269 ± 0.026 


2.2-2.4 


200 


1.251 ± 0.034 


2.4-2.6 


316 


1.350 ± 0.068 


Table 4: 


Average 


charge ratios 



disagree on the 4<r level. To reduce the discrepancy to about lcr we assume - in the spirit of 
the Particle Data Group p2| - that for all experiments an additional systematic error of ±0.015 
must be added, in form of a scale uncertainty common to all measurements of one experiment, 
independent of momentum. Clearly, this is a crude model! 

Including these errors results in the charge ratios displayed in figure [12|. Note that the values are 
quite close to those in figure ITUl while the error bars are enlarged. The corresponding numbers 
are listed in table f|. The momenta are calculated from the central values of the logarithmic 
bins. 

There is no final answer to question ii). While the measurements by Baxendale (and also others, 
with larger errors) indicate an increase of the charge ratio at momenta around 30 GeV, the data 
by Rastin do not support this hypothesis. For the moment the measurements are consistent 
with the simple hypothesis of a momentum independent charge ratio. 

We try to summarize the charge ratio measurements and their uncertainty (68% CL) with the 
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Figure 12: Summary of measured charge ratios and a parameterisation of the central value and 
the uncertainties. 



following formula: 



Rf, = 1.268 ± 



0.008 + 0.0002 



V 



GeV 



(17) 



in the momentum range 10 — 300 GeV. Figure shows the corresponding mean value and 
the error band. The central value is the mean of all measurements, taking into account the 
additional systematic error of 0.015. The momentum dependent error is estimated such that 
it is roughly of the same size as the uncertainties of the corresponding data points in figure [12 
and table [|. 

Theoretical models of atmospheric showers must be able to reproduce these data, the calculated 
charge ratios should fall into the band given in equation (O) and figure [Tj| Clearly, at high 
momenta precise data is still lacking. Momenta above a few hundred GeV are of particular 
interest, since a growing influence of kaons and a resulting increase of the charge ratio is 



predicted p3| . 



6. Summary and conclusions 

We have combined the published data on the vertical muon spectrum and charge ratio at sea- 
level. In this comparison we have found that the differential spectrum can be described using 
a simple formula. The shape of the momentum spectrum is well measured at momenta below 
100 GeV. Above 200 GeV only a few data points exist, therefore the uncertainty increases to 
17 % at 1 TeV. 

Several experiments measure the absolute normalization of the spectrum. Our combined result 
is compared to the calculation by Bugaev et al. At 10 GeV the measured flux is 11 % below 
the calculated one. 

The charge ratio is reported by many experiments. The combined result favors a momentum 
independent value of the charge ratio of 1.268. The error on the charge ratio increases rapidly 
above 200 GeV due to a lack of precise experimental data in that region. 
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APPENDIX A Measurements. 



The following three tables list the experiments/publications we considered, in chronological 
order, together with the most important parameters. 

Some entries are missing since the corresponding figures are not published. In particular the 
zenith angle regime is characterized frequently only by verbal expressions like 'near vertical'. 
We distinguish three magnet types: solid iron and conventional coil, air gap magnet with 
conventional coil and air gap magnet with superconducting coil. The period of data taking is 
indicated by year and month, e.g. 59/11 stands for November 1959. 

The following remarks refer to the experiment numbers in column 1 of the tables. The figures 
mentioned are those of the corresponding publication. 



1) Not used, since data set is a subset of 3) 

2) Spectrum data read off from figure 1; charge ratio taken from 3) 



6) Spectrum normalized to Rossi |28"1; data read off from figures 8 and 15. 

7) Spectrum normalized to Rossi. 

8) Spectrum normalized to Rossi; data read off from figures 1 and 3. 

9) The two values at p = 240 GeV are not statistically independent. We have calculated the 
arithmetic mean of the two figures and took the smaller of the two statistical errors as the 
uncertainty. 

10) We do not use these data, which were obtained in the equator region, where the geomagnetic 
cutoff is large (14 GeV) and may influence the muon flux. 

11) The spectrum data form a subset of those published in 13) jjTjJ]. Nevertheless we consider 
it separately, since the normalization procedures are slightly different. When calculating world 
averages we exclude these data. We do not use the charge ratio data, since they are included 
in the superset published in 13) |fl9[1 . 

13) Spectrum: no absolute flux determination, only normalization to previous measurements 
by other experiments. Charge ratio: A few measurements are given with slightly asymmetric 
errors; they have been 'symmetrized' by shifting the central value to the center of the error 
interval. 

15) Apparatus under concrete shelter of 868g/cm 2 . 

16) We use only the charge ratio value obtained at the town of Kiel[^0|; for the other mea- 
surements, made at muon momenta close to and below 10 GeV in the equator region, the 
geomagnetic cutoff is large (14 GeV) and may influence the muon flux. 

20) Resolution correction is based upon MDM = 100 GeV; if MDM of 350 GeV is used, spectrum 
is flatter and better consistent with 18) Points up to 100 GeV are considered reliable. 

21) Spectrum data are normalized to an integral intensity at 5 GeV. The spectrum provided is 
the theoretical spectrum that fits the data best. 

22) Above 50 GeV authors question results due to resolution. 

24, 25) Same magnet as 22)0. Charge ratio and its error read off from figure 2. 
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No author (s), name location coordinates alt. magnet detector(s) period MDM zenith spec, ratio 

reference /m / GeV 



1 


B.G. Owen and 
j. Li. Wilson, iy4y|ooj 


Manchester 
(Great Britain) 


53°N, 2°W 


50 


air 


Geiger, 
flash tubes 




30 


<8° 


no 


2 


D.E. Caro 
et ai, iyoU[izj 


Melbourne 
(Australia) 


38°S, 145°E 


50 


air 


counters 




« 50 




yes 


3 


B.G. Owen and 
J.Li. Wilson, iyoi|o4j 


Manchester 
(Great Britain) 


53°N, 2°W 


50 


air 


Geiger, 
flash tubes 




30 


<8° 


no 


4 


I. Filosofo 
et. ai, iyo4|ooj 


Agordo 
(italyj 


46°N, 12°E 


600 


iron 


counters 








no 


5 


B.G. Owen and 
Wilson, iyoo[ioj 


Manchester 
(Great Britain) 


53°N, 2°W 


50 


air 


counters 






< 10° 


yes 


6 


J. Pine et al, 
iyoy [i4j 


Cornell 


42°N, 76°W 


500 


air 


Geiger, 
cloud ch. 




175 




yes 


7 


J.E.R. Holmes 
et al, lybl[lbj 


Manchester 
(Great Britain) 


53°N, 2°W 


50 


air 


cloud ch., 
Geiger 


53-55 




< 10° 


yes 


8 


W. Pak et al, 
1961[15] 


Cornell 
(USA) 


42°N, 76°W 


500 


air 


Geiger, 
hodosc. 




175 




yes 


9 


P.J. Hayman and A.W. 
Wolfendale, 1962 [17, 10] 


Durham 
(Great Britain) 


54°N, 1°W 


70 


air 


Geiger, 
flash tubes 


59/11-60/03 


657 




yes 


10 


O.C. Allkofer 
et al, 1968 [37] 


near 
equator 


29°W, 0°N 
22°W, 1°S 





air 


spark ch. 








yes 



No author(s), name location coordinates alt. magnet detector(s) period MDM zenith spec, ratio 

reference /m / GeV 



11 


S.R. Baber 
et al, lyt>o|oo, oyj 




Nottingham 
(Great Britain) 


53°N, 1°W 


52 


iron 


Geiger 


64/06-65/05 


360 


yes 


yes 


12 


A.M. Aurela and A.W. 
Wollendale, iyo/[18J 




Durham 
(Great Britain) 


54°N, 1°W 


70 


air, 
iron 


Geiger, 
liasn 


64/06-65/01 




yes 


no 


13 


I.C. Appleton 
et al, iy i i[iyj 




Nottingham 
(Great Britain) 


53°N, 1°W 


52 


iron 


flash 
tubes 


64-68 


360 


yes 


yes 


14 


B.J. Bateman 
et al, iy / l|zlj 


AMH 


College Station 

/TTC A ^1 

(UoAJ 


31°N, 96°W 


80 


iron 


spark ch., 
scint . 






yes 


no 


15 


O.C. Allkofer 
et al, iy / l|zUJ 




Kiel 
(Germany) 


54°N, 11°E 


10 


iron 


spark ch., 
scint . 




> 1000 


yes 


no 


16 


O.C. Allkofer and 
W.-D. JJau, iyfz|4UJ 




Kiel 
(and equator) 


54°N, 11°E 


10 


iron 


spark ch., 
scint. 






no 


yes 


17 


B.C. Nandi 
et al, iy/z[zy, 22J 




Durgapur 
(India) 


24°N, 87°E 


70 


iron 


flash tubes, 
Geiger 


69/02-70/02 


985 


yes 


yes 


18 


C.A. Ayre 
et al, 1975 [23] 


MARS 


Durham 
(Great Britain) 


54°N, 1°W 


70 


iron 


scint., 
flash tubes 


72/05-73/01 


670 


yes 


no 


19 


J.M. Baxendale 
et al, 1975 [30] 




Durham 
(Great Britain) 


54°N, 1°W 


70 


iron 


scint., 
flash tubes 


72/02-72/12 




yes 


yes 


20 


P.J. Green 
et al, 1979[24] 


AMH 


Houston 
(USA) 


30°N, 95°W 


10 


iron 


spark ch., 
scint. 




345 < 9° 


yes 


no 



No author (s), name location coordinates alt. magnet detector (s) period MDM zenith spec, ratio 

reference /m /GeV 



21 B.C. Rastin Nottingham 53°N, 1°W 52 iron flash tubes, 74/09-78/05 3400 yes yes 
1984[25, 31] (Great Britain) scint. 

22 M.P. De Pascale MASS Prince Albert 53°N, 106° W 600 air MWPC, 89/08 118 yes yes 

et al, 1993[6] (Canada) (superc.) scint., TOF 

23 S. Tsuji Okayama 34°N, 134°E 5 iron drift, 92/09-97/12 270 0°-l° yes no 
et al, 1998[26] (Japan) scint. 

24 J. Kremer CAPRICE Lynn Lake 57°N, 101°W 360 air prop., drift, 94/07 175 0°-20° yes yes 
et al, 1999[27] (Canada) (superc.) TOF, RICH 

25 J. Kremer CAPRICE Fort Sumner 34°N, 104°W 1270 air drift, 97/04-97/05 175 0°-20° yes yes 
et al, 1999[27] (USA) (superc.) TOF, RICH 



APPENDIX B Muon flux data. 

The following lists contain all spectrum data 
for momenta above 10 GeV. The flux is given 
in (m~ 2 sr _1 s _1 GeV). Each line contains the 
momentum in GeV together with the pub- 
lished value and uncertainty of the flux. Al- 
titude corrections have been applied later and 
are not included in the figures listed here. 



J.E.R. Holmes et al, 1961 (111 



D.E. Caro et al, 1950 [jj 



12.2 
16.9 
17.8 
26.5 
51 



(6.46 ± O^-IO- 1 
(2.59 ± O^j-lO- 1 
(2.38 ± O^j-lO- 1 
(7.4 ± 0.9)-10- 2 
(1.3 ± 0.2)-10- 2 



11 


(7.88 ± 





32)-10~ 1 


13 


(5.90 ± 





30)-10 1 


16 


(3.68 ± 





lsj-io- 1 


19 


(2.33 ± 





12)- 10 1 


23 


(1.67 ± 





12)-10 _1 


28 


(1.04 ± 





07)-10- 1 


36 


(5.75 ± 





40)-10- 2 


49 


(2.02 ± 





20)-10~ 2 


67 


(8.65 ± 


1 


30)-10~ 3 


89 


(4.38 ± 





70)-10- 3 


134 


(1.42 ± 





26)-10" 3 


271 


(2.8 ± 





6)-10" 4 


1160 


(4.8 ± 


2 


3)-10~ 6 



B.G. Owen and J.G. Wilson, 1955 13 



P.J. Hay man and A.W. Wolfendale 
19620 



10.0 1.09 ± 0.03 
15.0 (4.36 ± 0.11)-10" 
20.0 (2.20 ± 0.06)-10" 



J.Pine et al, 1959 



11.6 


(8.26 ± 





99) 


10 


13.7 


(5.24 ± 





84) 


10 


16.3 


(3.05 ± 





49) 


10 


18.1 


(2.79 ± 





33) 


10 


19.8 


(1.91 ± 





31) 


10 


25.5 


(1.16 ± 





19) 


10 


29.8 


(7.68 ± 





92) 


10 


35.4 


(4.95 ± 





79) 


10 


46.0 


(1.99 ± 





32) 


10 


59.1 


(9.29 ± 


1 


49) 


10 


83.2 


(4.60 ± 





87) 


10 


125. 


(1.34 ± 





33) 


10 


et al, 


1961 |I§ 








12.9 


(6.14 ± 





68) 


10 


15.8 


(3.35 ± 





37) 


10 


19.8 


(1.89 ± 





21) 


10 


24.1 


(1.35 ± 





19) 


10 


26.4 


(1.03 ± 





11) 


10 


37.4 


(3.77 ± 





41) 


10 


58.9 


(1.29 ± 





14) 


10 



10.8 


(8 


51 


± 


0.26) 


■10- 


-1 


12.4 


(6 


14 


± 


0.17) 


■10- 


-1 


14.6 


(4 


35 


± 


0.11) 


•10- 


-1 


17.8 


(2 


52 


± 


0.07) 


•10- 


-1 


22.6 


(1 


39 


± 


0.04) 


■10- 


-1 


31.3 


(5 


85 


± 


0.15) 


•10- 


-2 


42.3 


(2 


88 


± 


0.11) 


•10- 


-2 


56.1 


(1 


22 


± 


0.05) 


•10- 


-2 


72.5 


(5 


75 


± 


0.46) 


•10- 


-3 


88.1 


(3 


27 


± 


0.29) 


•10- 


-3 


112 


(1 


36 


± 


0.14) 


•10- 


-3 


153 


(5 


18 


± 


0.93) 


•10- 


-4 


244 


(1 


14 


± 


0.22) 


•10- 


-4 


413 


(1 


98 


± 


0.99) 


■10- 


-5 


894 


(1 


84 


± 


1.01 


•10- 


-6 



A.M. Aurela and A.W. Wolfendale 



1967[|Tg 



15.1 (4.25 ± O.^-IO" 1 
41.5 (3.40 ± 0.44)-10~ 2 
82.1 (4.10 ± 0.35)-10-' 3 



25 



S.R. Baber et al, 19681551 



B.J. Bateman et al, 1971 |EI 



11.60 


(7.77 ± 0.26) 


io- 1 


10.0 


1.12 


15.22 


(4.22 ± 0.21) 


io- 1 


13.0 


(6.63 ± 


19.20 


(2.42 ± 0.12) 


io- 1 


15.0 


(4.56 ± 


24.00 


(1.39 ± 0.07) 


io- 1 


20.0 


(2.31 ± 


33.5 


(5.78 ± 0.35) 


io- 2 


27.0 


(1.08 ± 


50.0 


(1.90 ± 0.16) 


io- 2 


35.0 


(5.45 ± 


81.0 


(4.59 ± 0.60) 


IO" 3 


40.0 


(3.78 ± 


127.0 


(1.14 ± 0.18) 


10~ 3 


46.0 


(2.55 ± 


266.0 


(1.00 ± 0.24) 


10~ 4 


53.0 


(1.70 ± 


810.0 


(2.11 ± 0.55) 


10~ 6 







)-10 
)-10 
)-10 
)-10 
HO 
)-io- 2 
)-io- 2 
)-io- 2 



I.C. Appleton et al, 1971 IS 



O.C. Allkofer et al, 1971 



12.84 


(2.92 


± 


0.04) 


•10' 


_ 1 


11.4 


1.13 ± 0.01 


17.2 


(1.46 


± 


0.02) 


•10' 


_ 1 


14.8 


(6.04 ± O^-IO" 1 


24.3 


(5.78 


± 


0.13) 


■10' 


_2 


20.5 


(2.51 ± 0.03)-10^ 


33 


(2.54 


± 


0.06) 


•10' 


o 

— z 


31.4 


(8.01 ± 0.13)-10~ 2 


AO A 

4o.4 


( 1 nn 
(1.09 


± 


0.10 J 


i n- 
•10 


-2 


52.3 


(1.89 ± 0.05)-10~ 2 


A K P. 

45.0 


( 1 (Y7 

(l.UY 


i 

± 


n ^^^ 

u.uy J 


t n- 
•1U 


-2 


93.0 


(3.38 ± 0.14)-10- 3 


4o.z 


( d nn 
(y.UU 


± 


n qi "\ 
U.ol J 


t n- 
•1U 


-3 


175.0 


(5.19 ± 0.37)-10- 4 


51.0 


(6.95 


± 


0.67) 


■10' 


_3 


329.0 


(7.84 ± 1.12)-10- 5 


54.2 


(6.84 


± 


0.63) 


•10' 


-3 


642.0 


(6.40 ± 1.92)-10~ 6 


57.9 


(5.20 


± 


0.52) 


■10' 


-3 






62.0 


(4.79 


± 


0.47) 


■10' 


-3 


• B.C. Nandi and M.S. Sinha, 1972 


66.9 


(2.80 


± 


0.33) 


■10' 


-3 






72.6 


(2.76 


± 


0.31) 


■10' 


-3 


11.8 


(9.43 ± O.IS)-^" 1 


79.3 


(2.02 


± 


0.23) 


■10' 


-3 


14.0 


(6.38 ± o.uyio- 1 


87.5 


(1.57 


± 


0.17) 


■10' 


-3 


16.4 


(4.21 ± o.ogj-io- 1 


97.6 


(1.17 


± 


0.13) 


■10' 


-3 


19.7 


(2.72 ± o.oej-io- 1 


110.3 


(6.97 


± 


0.94) 


■10' 


-4 


24.2 


(1.41 ± 0.04)-10- 1 


129.9 


(3.00 


± 


0.54) 


■10' 


-4 


29.6 


(i.oi ± o.osj-io- 1 


149.4 


(2.93 


± 


0.29) 


•10" 


-4 


37.1 


(5.40 ± 0.17)-10~ 2 


181.5 


(1.34 


± 


0.24) 


•10' 


-4 


46.9 


(2.99 ± 0.13)-10- 2 


230.9 


(6.61 


± 


1.41) 


•10' 


-5 


60.0 


(1.45 ± 0.07)-10~ 2 


316.1 


(2.72 


± 


0.59) 


•10' 


-5 


84.0 


(5.58 ± 0.28)-10- 3 


491.5 


(7.67 


± 


1.69) 


■10' 


-6 


118 


(2.04 ± 0.17)-10~ 3 


1000.0 


(5.21 


± 


1.34) 


■10' 


-7 


167 


(6.09 ± 0.61)-10~ 4 














260 


(1.96 ± 0.25)-10~ 4 














467 


(2.69 ± 0.67)-10~ 5 














1109 


1.03 ± 0.36 -IO" 6 
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C.A. Ayre et al, 1975 M 



P.J. Green et al, 1979 



21 3 

Zj X . O 
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I Zj .UC/U 
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xu 
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xu 
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h 708 
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M.P. De Pascale et al, 1993 f 
Positive Muons: 
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(1.246 ± o.onj-io- 1 
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(3.09 ± 0.07)-10- 1 
(1.08 ± O.Osj-lO- 1 
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(7.1 ± 0.6) 

(3.0 ± 0.3) 
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10" 2 

io- 2 

10~ 3 
10~ 3 
10~ 3 



J. Kremer et al, 1999 
Positive Muons: 



1997 data 



12.42 
18.85 
26.68 
36.69 
51.47 
72.08 
100.96 



(4.14 ± O^-IO- 1 
(1.54 ± 0.04)-10- 1 
(6.4 ± 0.2) -IO" 2 
(2.8 ± 0.1) -IO" 2 
(10.2 ± 0.5) -IO" 3 
(4.2 ± 0.3) -IO" 3 
(1.5 ± 0.1) -10- 3 



Negative Muons: 



12.42 
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36.69 
51.47 
72.08 
100.96 



(3.20 ± 0.07)-10- 1 
(1.16 ± O.Osj-lO" 1 
(4.5 ± 0.2) -IO" 2 
(2.03 ± 0.08)-10- 2 
(7.7 ± 0.4) -10- 3 
(3.2 ± 0.2) -IO" 3 
(1.1 ± 0.1) -10" 3 



J. Kremer et al, 19990 1994 data 
Positive Muons: 
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APPENDIX C Charge ratio data. 

The following lists contain all charge ratio data 
for momenta above 10 GeV. Each line contains 
the momentum in GeV together with the pub- 
lished value and uncertainty of the charge ra- 
tio. 



P.J. Hayman and A.W. Wolfendale, 
1962 fill 



D.E. Caro et al, 1950(12 



35.0 1.6 ± 0.2 



B.G. Owen et al, 1951 134] 
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M.P. De Pascale et al, 1993 § 



19.89 
55.87 



1.292 ± 0.075 
1.409 ± 0.173 



J. Kremer et al, 1999 E7\ 1994 data 



12.42 
18.85 
26.68 
36.69 
51.47 
72.08 



1.257 ± 0.038 
1.269 ± 0.055 
1.372 ± 0.086 
1.466 ± 0.110 
1.384 ± 0.144 
1.212 ± 0.169 



100.96 1.235 ± 0.233 



J. Kremer et al, 1999071 1997 data 



12.42 1.291 ± 0.027 

18.85 1.335 ± 0.038 

26.68 1.427 ± 0.058 

36.69 1.383 ± 0.068 
51.47 1.320 ± 0.090 
72.08 1.337 ± 0.120 
100.96 1.337 ± 0.169 



APPENDIX D Correlation matrix for the fit of the spectrum shape: 





H 3 


H 2 


#i 


s 2 


#3 


1.000 


0.319 


-0.298 


-0.771 


H 2 


0.319 


1.000 


-0.105 


-0.558 


#1 


-0.298 


-0.105 


1.000 


0.672 


s 2 


-0.771 


-0.558 


0.672 


1.000 
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